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(57) Abstract: The described method and system are adapted to reduce the error between an ideally expected output signal and the 

actual output signal. The proposed adaptation algorithm is able to minimise, for instance in a system with a given transfer function, 
^ the error y-x between y=g(f(x)) and x, where g is an unknown and/or time-varing function and f the adaptive function for which the 

characteristic is changed to track g. The proposed adaptation algorithm updates not only the transfer function f at the current input 
^5 value x, but also the transfer function f at other points corresponding to different input values. One of the applications for such an 

algorithm is digital predistortion where a transmitter's non-linear characteristic needs to be linearised, in an adaptive manner, since 
^ the characteristic exhibits slow changes with temperature, bias, ageing or the like. 
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Method and system for compensating non-linearities and time- 
varying changes of a transfer function acting on an input 
5 signal 

Description 
10 Field of the Invention 

The invention is directed to a method for compensating non- 
linearities and time-varying changes of a transfer function 
acting on an input signal. This input signal is first subjected 

15 to a first, adaptive transfer function and thereafter to a 

second, unknown and varying transfer function to generate the 
actual output signal. When the second transfer function varies, 
the first transfer function is updated for compensating these 
changes. The described method and system are therefore adapted 

20 to reduce the errors between an ideally expected output signal 
and the actually generated output signal. 

Background of the Invention 

25 

For example, in the field of mobile communications, there is 
sometimes the need to linearise a transmitter's non-linear 
characteristic, in an adaptive manner, when the characteristic 
exhibits slow changes caused by temperature, bias, ageing or 

30 the like. One of the existing compensating possibilities is the 
digital predistortion of the input signal before applying same 
to the transmitters power amplifier. Thereby, the error 
between an ideally expected output signal and the output signal 
actually generated in response to a current input signal can be 

35 minimised. This need for compensation not only occurs in the 

field of mobile communications but also in other fields where a 
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transfer function varying in an unknown manner is to be 
compensated by adapting an adaptive transfer- function. 

When the characteristic is changing, it takes some time until 
5 the system has been adapted to the new situation. It is 
desirable to perform this adaption as swiftly as possible. 
However, normally, a great number of iterations is necessary 
until the system has actually been adapted to the new 
characteristic. 

10 

Summary of the Invention 

It is therefore an object of the present invention to provide a 
15 method and system which are able to swiftly adapt to a changing 
transfer function and thereby quickly minimise an error between 
the actual output signal caused by a current input signal, and 
an ideally expected output signal. 

20 The invention provides a method for compensating deviations of 
an unknown transfer function from an expected transfer 
function, and/or for compensating time-varying changes of a 
transfer function acting on an input signal, for minimising 
errors of an output signal generated in dependence on the input 

25 signal, which input signal is subjected to a first, adaptive 

transfer function and to a second, varying transfer function to 
generate the output signal, the first transfer function being 
updated for compensating deviations of the second transfer 
function from an expected transfer function and/or for 

30 compensating changes of the second transfer function, wherein, 
when updating one point of the first transfer function for a 
current input signal value, the first transfer function is also 
updated for at least one other point corresponding to a 
different input signal value. 

35 

The invention furthermore provides a system for compensating an 
unknown and/or varying transfer function acting on an input 
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signal, for minimising errors of an output signal generated in 
dependence on the input signal, which input signal is applied 
via a first, adaptive compensating means having a first, 
adaptive transfer function, to a second means having a second, 
5 unknown and/or varying transfer function to generate the actual 
output signal, the first transfer function being updated for 
correcting deviations of the second transfer function from an 
expected (ideal, wanted) transfer function by means of a 
processing means, wherein the processing means is adapted to 
10 update, when updating one point of the first transfer function 
for a current input signal value, the first transfer function 
also at at least one other point corresponding to a different 
input signal value. 

15 Due to this multi-point-correction of the first transfer 

function, the described method and system are able to quickly 
adapt the first transfer function to changes of the second 
transfer function and/or to quickly adapt the overall transfer 
function to the desired one. Thus, the apparatus swiftly 

20 converges to the new (or desired) condition, with a drastic 
reduction of the necessary adaption time. The occurrence of 
errors between an output signal ideally to be expected for a 
given input signal, and the actually generated output signal is 
therefore limited to a very short time interval after a 

25 variation (or the first use) of the second transfer function. 

Furthermore, the error deviations within this time interval are 
reduced to smaller values. 

30 Brief Description of the Drawings 

The present invention will be more readily understood upon 
referring to the following description of preferred embodiments 
when read in conjunction with the accompanying drawings. In the 
35 drawings, 



Fig. 1 shows a basic configuration of a functional system, 
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Fig. 2 shows a chart for explaining the algorithm used in the 
present invention, 

Fig. 3 shows the system of Fig. 1 in a practical implementation, 

Fig. 4 shows an iterative updating of two values of a transfer 
function, 

Fig. 5 illustrates the updating of a plurality of values of a 
transfer function, 

Fig. 6 shows an embodiment of the invention implemented as an 
adaptive predistorter and amplifier, 

Fig. 7 shows an example of the magnitude and phase response of 
function g, 

Fig. 8 illustrates an example of the ideally expected and 
actually generated output signals for a non-linearised function 

7* 

Figs. 9 and 10 show simulation results for the presented 
algorithm and known algorithms, 

Fig. 11 illustrates F values calculated at the end of the 
simulation, 

Figs. 12 and 13 show the in-phase and quadrature-phase 
components of input and output signals, 

Figs. 14 and 15 illustrate the different convergence behaviour 
of the presented and known methods for different variations of 
the transfer function g, and 

Figs. 16 and 17 show the results of three different updating 
procedures usinq the linear and LMS methods, resp. 
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Detailed Description of Preferred Embodiments of the Invention 

Generally, the described adaptation algorithm is a method of 
minimising the error |y-x| between y=g(f(x)) and x, where g is 
an unknown function and f the adaptive function for which the 
characteristic is changing to track g. Ideally, 
f=g~ 2 . Unlike current algorithms, it updates not only the 
transfer function f at the current input value x, but also the 
transfer function f at other points (corresponding to different 
input values x updat e...) . One of the applications for such an 
algorithm is digital predistortion where a transmitter' s non- 
linear characteristic needs to be linearised, in an adaptive 
manner since the characteristic exhibits slow changes with 
temperature, bias, ageing.... 

First, the new . algorithm used in the method and system will be 
mathematically described. 

Let us consider an unknown function g defined as: 
g: D g c C -> C 

x g(x) 

C designates the complex domain, and can be extended to an n- 
dimension domain; D g : definition domain of g. 

The adaptive function f is defined as: 
f : Df c C C Df : definition domain of f 

x -» f(x) C: complex numbers 

The functions (transfer functions) f and g are performed in the 
order shown in Fig. 1. Fig. 1 illustrates a functional system 
having a compensating means 1 to which an input signal x is 
supplied at a first input terminal. The compensating means 1 
generates an output signal z in accordance with its transfer 
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function f (correspondence between x and z) .and the value of x. 
The signal z is input into a second means 2 which outputs the 
output signal y and has a transfer function g which may change 
over time. The output signal y is fed back to a second input 
5 terminal of the compensating means 1. As the overall transfer 
function is to be maintained constant, any change of the 
transfer function g or deviation from the ideal transfer 
function will be compensated by appropriately adapting the 
transfer function f so as to maintain the desired relationship 

10 between x and y. The means 1 may for example be an adaptive 

predistorter whereas means 2 may be a power amplifier. When the 
desired relationship between x and y is equity (amplification 
gain factor = 1) for example, any deviation of the gain of the 
amplifier 2 from the value 1 will be compensated by setting the 

15 amplification characteristic of means 1 to the reciprocal 
value . 

The algorithm attempts, in the present example with unity gain 
(gain factor = 1) , to minimise the error |y-x| between the 
20 output y and the input x, in an iterative process (see Fig. 2) . 

Fig. 2 is a diagram showing the characteristic (transfer 
function) g and the desired linearised overall response 
characteristic 3 (x=y) . The horizontal axis represents the 
25 absolute value of the output z of means 1 whereas the vertical 
axis represents the final output y (absolute value) . 

In the following, f k {x) denotes f(x) at iteration k, where an 
iteration is the process of updating simultaneously one or more 
30 values defining the function f. It will also be written as: 
f k (x)=x.F k , x (with F k , x eC) . 

At iteration 0 (k=0) , the transfer function f is set to 1: 
VxeC, F 0 ,x=l+j-0 
35 (i.e., VxeC, f 0 (x)=x). 
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y = g(z) = g(f k (x)) = g(x.F k , x ) = x.F k , x .G k 



(1) 



with x . F k , x =fk (*) and G 



y = g(xJ ? ktX ) 
z x.F kx 



Ideally, we would like: 



V x€C, y = x 



substituting in (1), 
o x.F k , z .G k , x = * 



(2) 




(defined as such) 



(3) 



F op t would be the optimum value to use for F k , x to satisfy (2) . 
Unlike previously published algorithms (linear method, secant 
method, Least Mean Square (LMS) etc., which have been used more 
specifically for predistortion systems), it is not 
theoretically correct to update the new value F k+ i, x with F opt 
(even though adaptation factors are used to avoid 
instability). If F k+ i, x =F opt would be given, (2) would become: 
x.F k+ i fX .G k , x = x which is the desired result. 

However, this does not take into account the fact that changing 
F k+ i, x would mean x.F k+ i, x * x.F k , x , yielding a different input to 
g (see Fig. 2) . 

Let's define yi and y2 as: 

yi = y =g(fk(x)) = g(x.F k/X ) = x.F k , x .G k , x 

y 2 = g(fk+i(x)) - g(x.F k+ i, x ) = x.F k+ i fX .G k+1/X (with F k+ i, x =F opt here) 
therefore y 2 * x because G k , x * G k +i, x 

What really ought to be updated is the value of f corresponding 
to a different input x upda te, while F k+ i, x remains unchanged 
(F k+ i, x =F k(X ) . x U p dat e is defined as: 
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O Xupdate- F k+ i, xupdate 88 X.F k , x with F k+ i, xupdate ~ F opt 



<=> X update = — ^ J 



substituting for F opt from (3) yields, 

Xupdate = X.F k , x .G k ,x (5) 

comparing (5) with (1), 

<=> Xupdate = y (6) 

Therefore, now (iteration k+1) , x up date yields the same input to 
g as was x at the previous iteration (iteration k) : 

A+l ( X update ) = A ( X ) 

Z 

gCA + l(X update ))=£(A00) 

* y ' 

assuming g has not varied from 
iteration k to k+l 

<=> G k +i, xupdate = G k , x 

<=> yupdate = g(^k+l (Xupdate) ) = 9 (Xupdate ■ F k+ i, xupdate) = ^update • F opt ) 

<=> yupdate = g(x.F k , x ) from (4) 

(7) 

o yupdate = y W ' 

Recalling (6) and having (7), the result is 

( 8 ) 

yupdate ~ Xupdate 

which corresponds to a linearised system at point x upd ate 
(remember the aim was to have y=x) . 



In this form however, the algorithm may present the danger of 
not updating f completely since the current jf k , x is pointing at 
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another value f k+ i, xupdate to be updated. And once the latter has 
been updated/ it will not point to another fjc+i,x3 ( x 3 *x U pdate t 
determined in the same manner as x upd ate was obtained from x) but 
merely to itself. 

Bearing this in mind, the problem is preferably overcome by 
updating the current value f k+ i, x for the current input x too. 
The following overall solution is suggested: 



F k +i, xupdate = /*.F op t or similar (a) 
F k +i,x = chosen algorithm (linear, LMS, secant etc..) (b) 
where jueK (real numbers) and is constant (^adaptation f actor u ) . 

It is important to note that the key point here is to update 
not only the F k/X value corresponding to the current input x 
(i.e. F k +i,* - this is done in (b) ) , but also F k , XU p da te 
corresponding to a different value x U p da te defined by (4) (i.e. 
F k+ i, xupdate - this is done in (a)). Hence, a different algorithm 
could be used to update F k , x (F k+ i fX ) , like LMS, the linear 
method, Recursive Least Squares (RLS) etc... In other words, the 
important part is (a), whereas any algorithm can be used for 
(b) ; the adjunction of (a) to (b) improves (b) . 

(X F + J3 F 

One could choose for (b) : F k+ i, x = ' k,x - — with (a,fi)eR 2 . 

Note that, for example, if a=0=fJ=l and assuming the algorithm 
converges, the following holds true: 

F k ., L -» Fop, 



u,a,P can be chosen as a function of the input x, and therefore 
be different for different x inputs. 
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In the foregoing explanation, the equation x=y was assumed as 
desired relationship between x and y. In practise, however, the 
present algorithm can also be applied, with appropriate 
changes, to a more complicated dependency of y on x, for 
instance: y=A,.x (XeC, the previous description used . 

It can furthermore be envisaged to update adjacent indices to 
Fkn,xu P date as well (corresponding to other inputs), with 
appropriate weighting. This is because g can be assumed to have 
a continuous derivative (g is „smooth") ; doing this results in 
an even faster convergence (a particular illustration is 
discussed in the next section below) . 

One possible variation consists in updating f for all values 
between x and x upd ate, by interpolation. A possible way of doing 
this is to use linear interpolation. Assuming x < x up date/ 

r k+l,xtnlerpolnle " r k+!.x T _ V A interpolate > 

X update K 

withX huerpola« C e [ X ' X updnte] 

This can increase convergence speed quite considerably. 

In the following, a practical implementation will be described 
where D f is finite. 

In practise, the above described algorithm may be implemented 
in a computer program or in a DSP (digital signal processor), 
requiring digitisation. It is then suitable to limit the size 
of D f and define f as a list of complex values in a table which 
are constantly updated. The table can be multi-dimensional and 
the addressing customised to the application (e.g.: one could 
use the magnitude of the input value to address the table, or 
the square magnitude/power etc..) . 

The following example is an illustration of this particular 
case with a one-dimensional table of size m, addressed by the 
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magnitude of the input x. Since the table is finite, each input 
value x is associated with an index to the table 
(quantisation), fetching the corresponding value F ifk in the 
table. For example: f k (x) = x.F i#k ; f k (x) can also be produced 
by interpolation. 

Figure 3 shows a practical implementation of this system 
similar to same of Fig. 1 wherein the compensating means 1 is 
implemented as a table memory (here: one-dimensional look-up 
table) 4 which is addressed by the actual value of the input x. 
The values memorised in the table memory 4 are used for 
multiplication with the value of x for generating z, and are 
iteratively updated. 

Figure 4 illustrates the update process and shows, in the left 
part, the table contents (values Fi, k to F m , k ) at iterative step 
k, whereas the right part of Fig. 4 illustrates the table 
contents (values F 1>k+ i to F Bfk+ i) at the next iterative step k+1. 
The arrows between the left and right parts of Fig. 4 indicate 
the double updating not only of F i#k+ i but also of Fj^+i, 
corresponding to x and x upd atef resp. 

Applying the interpolation discussed above results in the 
diagram shown in Fig. 5. This can increase convergence speed 
considerably, however, it is computationally more complex and 
resource consuming. Fig. 5 illustrates the structure wherein 
more than two F values are updated based on the present value 
of x. Similar to Fig. 4, Fig. 5 shows, in the left part, the 
table contents (values F lfk to F m#k ) at iterative step k, whereas 
the right part of Fig. 5 illustrates the table contents (values 
Fi, k+ i to F m , k+ i) after the next iterative step k+1. Here, in 
addition to the updating not only of F i/k+ i but also of Fj, kn all 
values of F assigned to values of x lying between x and x upda te 
are updated (Fi-i, k+ i ... F j+x ,k+i) , for example by interpolation. 
The plurality of arrows between the left and right half parts 
of Fig. 5 indicate this multi-updating based on one actual x 
value. - . 
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The algorithm described above has to be adapted to this 
application, but basically remains very similar. We can choose 
for instance to address the F table by the magnitude of the 
5 input x (affects look-up table resolution) . 



It is assumed that the input x has a known limited magnitude 
range: |x|e[|x| min ; |x| max ]. The F table size is set accordingly 
to cover all the input values of x within this range. The total 

10 range of magnitude of input x is divided into subranges of 

possibly, but not necessarily equal size (it can for instance 
be a function of the input signal's statistical distribution 
. . . ) , the number of subranges corresponding to the number m of 
values of F contained in the table. Therefore, each actual 

15 value of x addresses that table value F which is attributed to 
that subrange to which the present x belongs, for instance Fj.. 
Mathematically: 

for | x ) g [ax, a i+ i ] , x ^ index i in the F table corresponding to ¥ ik 

associated to 

with aiGR + , and ie[l->m] 
20 Hence, f k (x) = x. Fi, k (interpolation may be used here to 
determine f k (x) ) . 



The index i can be determined as: 

12 N 



i = floor 



(m-l).|x| 2 



+ 1 



(9) 



25 Equation (9) makes sure that the indices i range from 1 for 
lx|=0 to i=m for |x|=|x| roa x- 



Repeating the algorithm as it has been described in part 1: 
y = g(fvAx)) = g(x.Fi, k ) = x. Fi, k .G if k,x 

y g(x.F itk ) 



30 with x.F irk = f k (x) and G ikx =- = 



x.E 
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As this practical implementation of the system has a finite D f , 
the additional index "i" is added which refers to the value F 
in the look-up table corresponding to cell "i" (see Fig. 4). 
"k" is the iteration number, whereas "x" is used as a referring 
index for the gain g at point x.F i/)c which is the input to the 
transfer function g. 

x F 

So x update - — — and the associated index in the F table for 
Xupdate is j, with j*i. 

It must be ensured that |x updat el is within the range covered by 
the F table (i.e., Ixlmin ^ lx upd atel ^ lx| max <=> j€[l-» 
m]) . 

If lx update l < lx| min then j is preferably set to j=l 
and if |x| max < |x upda tel then j is preferably set to j=m. 
Alternatively, the F table could simply not be updated at index 
j for j£ [l-»m] . 

Finally, the following settings are selected, similar to the 
above discussed case: 

Fj,k+i^.F op t (i) 
Fi,jc+i= chosen algorithm (linear, LMS, secant etc..) (ii) 

where /ieR and is constant. 

If i=j, one may choose either (i) or (ii) for the update. 

This algorithm has been simulated in a system commonly referred 
to in the literature as predistortion . 

Fig. 6 shows a predistortion block diagram. The input x is 
complex and has been represented as: 
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x(t)=sin(co.t)+j . sin(co.t) 



ca=2x7ix8.7 51e3 rad.s' 1 , 



|x(t)| = V2x|sin(6).t)j 

Zx(t)-fW 
4 



Sampling time=iteration f requency=200e3 Hz. 
F table size is m=128. 

(The expression "e3" means a multiplication by a factor 10 3 ) . 
In Fig. 6 and the following figures, the following applies: 

i in =Re(x); q in =Im(x) [input - x] 

i pd =Re(z); q pd =Im(z) [predistorter out - fix)] 

ife=Re(y); q fe =Im(y) [„PA out" - g(f(x))] 

(PA = power amplifier) . 

The adaptive predistorter amplifier shown in Fig. 6 comprises a 
multiplier 5 to which the input signal x and the output of a 
look-up table (memory) 6 is applied. The contents of look-up 
table 6 correspond to the table shown in Fig. 4 or 5. Thus, the 
output signal of look-up table 6 is the F ifk value selected 
according to the present amplitude of input signal x. The 
multiplier 5 generates, as its output, the signal z which is 
applied to means 2, similar as in Figs. 1 and 3. Here, the 
means 2 is an amplifier having a time-varying gain 
characteristic. The output signal y of means 2 and the input 
signal x are supplied to a processing means 7 which stores and 
processes the above described adaptation algorithm. Updates for 
the indices i, j, etc., i.e. F values F i/k+ i, F jfk +i, etc... for 
the present iteration step "k+1" are calculated within 
processing means 7 based on the algorithm and are then supplied 
from processing means 7 to the look-up table 6 for storage 
therein, replacing the previous respective values (see Fig. 4 
or 5) . 

The complex input x is not only supplied to processing means 7 
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but also to a squaring means 8 which generates the square value 
of the absolute value of input x. The output of squaring means 
8 is connected to an input of look-up table 6. The present 
output value of squaring means 8 is used as address (index i) 
for selecting the F ifk value to be updated, so that same is 
selected depending on the square of the absolute value of x. 
This selected F value and the at least one further F value (see 
Fig. 4 or 5) are then updated according to the calculation 
results of processing means 7. 

The characteristic of g (means 2) is shown in Fig. 7. The upper 
diagram of Fig. 7 represents the magnitude response of g 
(horizontal axis: absolute value of input signal z; vertical 
axis: output y= g(|z|)) whereas the lower diagram illustrates 
the phase response of g (horizontal axis: absolute value of 
input signal z; vertical axis: phase of g(|z|) in radians). 
These curves correspond to typical measured data for a given RF 
power amplifier (data has been normalised) . Typically, such an 
amplifier will produce essentially AM to AM and AM to PM 
distortion. 

For reference, Fig. 8 shows the input and output signals for a 
non-linearised function g (i.e., VxeC, f k (x)=x <=> V(i,k)e[l-> 
m]xN, Fi, k (x)=l+j .0 => y=g(f (x) )=g(x) ) with predistortion 
being switched off. The nipper diagram of Fig. 8 illustrates the 
time-behaviour of the real parts of input x (i in = Re(x)) and 
output y (i fe = Re(y)) with regard to the number of iterations 
(horizontal axis). The lower diagram of Fig. 8 illustrates the 
time-behaviour of the imaginary parts of input x (q in = Im(x)) 
and output y (q fe = Im(y)) with regard to the number of 
iterations (horizontal axis) . Here, we are aiming at a gain 
factor of 1 so the input signals are similar to the ideally 
expected output signals. 

Let us first consider a case where the function g is unknown 
but time-invariant. Its transfer function is shown in Fig. 7. 



r f 
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for the presented algorithm. The upper diagram shows the 
imaginary part of input x over the number of iterations whereas 
the lower diagram shows the imaginary parts of signals z and y. 

5 It can be seen from the simulations that the presented 

algorithm presents faster convergence than the linear method or 
the LMS method. 

Next, let us consider a case wherein g increases monotonically 
10 with time. 

The new unknown transfer function g lt2 is defined as a ramp of 
slope 1 multiplied by the previous transfer function g: 
gi.2 (t, x)=g(txx) where g is the same function as defined in Fig. 
15 7. 

Again, Fig. 14 presents the results for this time varying 
function g (ramp) for all three methods, with the inscriptions 
of the horizontal and vertical axes similar to same of Figs. 9, 
20 10. 

Here again, the proposed algorithm reveals to converge faster 
and hence to be more robust against slow variations of the 
function g (gain) . 

25 

In the next considered case, the transfer function g is a time 
varying function g in form of a step. This new transfer 
function g U3 , is defined as the previous function g multiplied 
by a step function: 

30 

gi.3(t,x) = g(lxx) for 0<t<0.01 (t=0.01 <=> iteration 2000) 
gi. 3 (t,x) = g(1.2xx) for 0.01<t 



35 Fig. 15 presents the results for all three methods. 
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This case attempts to test the response of the algorithm for 
sudden variations of g. The proposed method again converges the 
quickest. Note that the LMS method is quite slow to readapt. 

A comparison for g as a step (same as above) is done for the 
linear method and the LMS method for the three following 
algorithms (Figs. 16 and 17): 

1) Updating current x address only (index i); 

2) Updating current x (index i) and x upda te (index j) 
addresses; 

3) Updating current x (index i) , x upda te (index j) and all 
xinterpoiate addresses between x and x upda te (indices w with we 
[ j->i] if j<i for example) . 

The algorithm 1) is the known one whereas algorithms 2) and 3) 
are in accordance with the present invention. 

Fig. 16 shows the three updating procedures using the linear 
method, whereas Fig. 17 shows the three updating procedures 
using the LMS method. The upper diagram illustrates algorithm 
1), the middle diagram algorithm 2), and the lower diagram 
algorithm 3) . 

It is obvious that algorithm 2) has a better convergence speed 
than algorithm 1) , and that algorithm 3) yet improves this 
speed dramatically. Interpolation used in algorithm 3) relies 
on the assumption of a „well-behaved" response for g, which 
becomes even more plausible as the algorithm converges. In 
general, for a power amplifier, this is a fair assumption. 

The algorithm described here thus provides, compared to other 
algorithms used for predistortion systems such as the linear 
method or the LMS method, very fast convergence, and is quite 
robust against fluctuations of g in time. Besides, it requires 
few operations (i.e. is computationally simple), and can be 
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used in conjunction with other algorithms as an enhancement. 

The use of spectrally more efficient modulation schemes (higher 
throughput/occupied bandwidth ratio) to accommodate higher data 
5 rates in communications often introduces both phase and 

magnitude variations of the signal, thus requiring stringent 
linearity of the transceiver to preserve signal integrity. 
Transmitter linearisation, in particular, is then often 
necessary to provide good linearity while offering reasonable 
10 power efficiency. 

Satellite and mobile applications are two potential candidates 
for linearisation. The proposed algorithm can, for instance, be 
used within digital predistortion systems as a linearisation 

15 means. The present invention can therefore be embodied in a 
stationary or mobile device for mobile communications, in 
particular in a RF transmitter thereof, in such standards as 
GSM EDGE (Global System for Mobile Communications Enhanced Data 
. Rates for GSM Evolution) , TETRA (Terrestrial Trunked Radio) , 

20 CDMA (Code Division Multiple Access), W-CDMA (Wideband-CDMA), 
NADC (North American Digital Cellular), PDC (Personal Digital 
Cellular) . . . 



An embodiment is therefore a stationary or mobile apparatus for 
25 mobile communications containing a device as shown in Fig. 3 or 
6, preferably as defined in the system claim, and/or operating 
as defined in the independent method claims. 
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1. Method for compensating deviations of an unknown transfer 
function from an expected transfer function, and/or for 
compensating time-varying changes of a transfer function acting 
on an input signal, for minimising errors of an output signal 
generated in dependence on the input signal, which input signal 
is subjected to a first, adaptive transfer function and to a 
second, unknown and/or varying transfer function to generate 
the output signal, the first transfer function being updated 
for compensating deviations or changes of the second transfer 
function, 

wherein, when updating one point of the first transfer function 
for a current input signal value, the first transfer function 
is also updated for at least one other point corresponding to a 
different input signal value. 

2. Method according to claim 1, wherein the at least one other 
point of the first transfer function is calculated in 
dependence on the current input signal value and the value of 
the first transfer function at this input signal value. 

3. Method according to claim 1 or 2, wherein the at least one 
other point of the first transfer function is calculated in 
additional dependence on the respective value of the second 
transfer function. 

4. Method according to claim 1, 2, or 3, wherein the at least 
one other point of the first transfer function is corrected in 
dependence on the value of the second transfer function for the 
current input signal value. 
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5. Method according to claim 4, wherein the at least one other 
point of the first transfer function is corrected by 
multiplying the value of the second transfer function for the 
current input signal value, with a factor. 

6. Method according to any one of the preceding claims, wherein 
values of the first transfer function lying adjacent to the one 
point of the first transfer function and/or the at least one 
other point of the first transfer function are also corrected 
with appropriate weighting. 

7. Method according to claim 6, wherein values of the first 
transfer function lying between the one and the at least one 
other point of the first transfer function are corrected by 
interpolating. 

8. Method according to any one of the preceding claims, wherein 
the values of the first transfer function are provided in form 
of a table being addressed by the input signal values. 

9. Method according to any one of the preceding claims, wherein 
the updating of the first transfer functio n is done in an 
iterative manner. 

10. System for compensating deviations of an unknown transfer 
function from an expected transfer function, and/or for 
compensating time-varying changes of a transfer function acting 
on an input signal, for minimising errors of an output signal 
generated in dependence on the input signal, which input signal 
is applied via a first, adaptive compensating means having a 
first, adaptive transfer function, to a second means having a 



WO 01/56146 PCT/EP00/00579 

22 

second, unknown and/or varying transfer function to generate 
the actual output signal, the first transfer function being 
updated for compensating deviations and/or changes of the 
second transfer function by means of a processing means, 
5 wherein the processing means is adapted to update, when 

updating one point of the first transfer function for a' current 
input signal value, the first transfer function also at at 
least one other point corresponding to a different input signal 
value. 

0 

11. Use of the method according to anyone of claims 1 to 9, or 
of the system according to claim 10, to perform linearisation 
of a signal path. 
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Figure 1: The functional system 




Figure 2: Illustration of the algorithm 
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Figure 3: Practical system 




Figure 4: F table 
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Figure 5: Variation for F table update 
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Figure 6: Predistortion block diagram 
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Figure 7: Magnitude and phase response of g 
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Figure 8: Non-linearised function g - (predistortion off) 
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Figure 9: 20*logl0(|y-x|) for the presented algorithm 
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Figure 10: 20*logl0(|y-x|) for linear and LMS methods 
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Figure 12: In-Phase components through the chain (input, predistorter out, PA out) for the 
presented algorithm 
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Figure 13: Quadrature components through the chain (input, predistorter out, PA out) 
the presented algorithm 
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Figure 15: time varying function g (step) 
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Figure 16: Three updating procedures for the linear method 
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Figure 17: Three updating pr cedures for the LMS method 
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